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ABSTRACT 

The influence of a dust grain mixture consisting of spherical dust grains with 
different radii and/or chemical composition on the resulting temperature struc- 
ture and spectral energy distribution of a circumstellar shell is investigated. The 
comparison with the results based on an approximation of dust grain parameters 
representing the mean optical properties of the corresponding dust grain mixture 
reveal that (1) the temperature dispersion of a real dust grain mixture decreases 
substantially with increasing optical depth, converging towards the temperature 
distribution resulting from the approximation of mean dust grain parameters, and 
(2) the resulting spectral energy distributions do not differ by more than 10% if 
> 2^ grain sizes are considered which justifies the mean parameter approxima- 
tion and the many results obtained under its assumption so far. Nevertheless, 
the dust grain temperature dispersion at the inner boundary of a dust shell may 
amount to ^100 K and has therefore to be considered in the correct simulation 
of, e.g., chemical networks. 

In order to study the additional influence of geometrical effects, a two- 
dimensional configuration - the HH 30 circumstellar disk - was considered, using 
model parameters from Cotera et al. (2001) and Wood et al. (2002). A drastic 
inversion of the large to small grain temperature distribution was found within 
the inner ~ lAU of the disk. 

Subject headings: radiative transfer — scattering — methods: numerical — cir- 
cumstellar matter 



1. Introduction 

The simulation of spectral energy distributions (SEDs), images, and polarization maps 
of young stellar objects has become a profound basis for the analysis and interpretation of 
observing results. Many techniques and approximations for the solution of the radiative 
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transfer (RT) problem in different model geometries (1D-3D), considering more and more 
special physical processes, such as the stochastic and photo-electric heating of small grains 
(see, e.g, Draine & Li 2001; Bakes & Tielens 1994, Siebenmorgen, Kriigel, & Mathis 1992), 
scattering by spheroidal grains (see, e.g.. Wolf, Voshchinnikov, & Henning 2002; Gledhill & 
McCall 2000), or the coupling of line and continuum RT (see, e.g., Rybicki & Hummer 1992), 
have been developed. The simulation of the temperature structure in simple-structured cir- 
cumstellar shells or disks (see, e.g., Malbet, Lachaume, & Monin 2001; Chiang et al. 2001), 
the estimation of the properties (luminosity, temperature, mass) of heavily embedded stars 
(see, e.g., Kraemer et al. 2001), or the determination of the inchnation of a circumstellar 
disk (see, e.g., Chiang & Goldreich 1999; Men'shchikov, Henning, & Fischer 1999; Wood 
et al. 1998) represented modest, first attempts of the application of the existing sophisti- 
cated numerical techniques. More recent efforts are directed to derive the dust grain size 
distribution in a circumstellar disk from its SED (Wood et al. 2002; D'Alessio, Calvet & 
Hartman 2001). Thus, it is clear that beside strong observational constrains, the model 
parameters and considered physical processes have to be questioned in depth in order to 
derive such detailed conclusions. However, looking behind the scenes, many of the RT mod- 
els are based on simplifying assumptions of very basic processes such as isotropic instead 
of anisotropic scattering, mean dust parameters representing dust grain ensembles (different 
radii and chemical compositions), or the flux- limited diffusion approximation - approxima- 
tions which are well-suited for handling the energy transfer in hydrodynamic simulations or 
a rough data analysis but which may not necessarily guarantee the desired accuracy for a 
detailed SED and image/polarization data analysis. 

In the advent of (space) observatories such as SIRTF^, which will be able to obtain 
SEDs of evolved debris disks around young stars, there exists a strong need for adequate 
numerical RT techniques in order to allow to trace dust grain growth (Meyer et al. 2001) and 
the influence of other physical effects and processes such as the Poynting-Robertson effect 
(see, e.g., Srikanth 1999) and dust settling (DubruUe, MorfiU, & Sterzik 1995, Miyake & 
Nakagawa 1995) on the dust grain size distribution in these disks. Therefore, we clearly have 
to understand which influence the different approximations (as long as they are required) in 
the RT simulations may have on the resulting observables. Based on two different grain size 
distributions consisting of astronomical silicate, the differences in the resulting dust grain 
temperature distributions and the resulting SEDs between simulations of the RT in (A) 
a "real" dust grain mixture and (B) under the assumption (approximation) of mean dust 
grain parameters will be discussed. This investigation is therefore mainly focused on two 
questions: (1) Of what order of magnitude are the differences (?) and (2) How many grain 
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sizes have to be considered to represent the properties of a real grain size distribution? The 
spatial temperature distribution of the considered dust configurations is calculated on the 
basis of local thermal equilibrium. Stochastic heating processes which are expected in case 
of very small grains consisting of tens to hundreds of atoms (see, e.g., Draine & Li 2000 and 
references therein) are not subject of this investigation. 

In Sect. 2, the RT and the dust grain model are briefly introduced. In Sect. 3, the 
definition of the mean dust grain parameters is given and the expected deviations of the RT 
results - once based on the mean dust parameters and once on a real grain size distribution 
- arc outlined. In Sect. 4, the RT in a spherical shell with variable optical depth and density 
distribution (see Sect. 4.1) is considered, while the temperature structure in a model of the 
HH30 circumstellar disk is investigated in Sect. 4.2. The SED resulting from a models with 
a dust grain mixture is compared to the mean particle approximation in Sect. 5. 

2. Radiative transfer and Dust model 

The RT simulations presented in this article have been performed with the three- 
dimensional continuum RT code MC3D which has been described by Wolf & Henning 2000 
(see also Wolf, Henning, & Stecklum 1999; Wolf 2002). It is based on the Monte-Carlo 
method and solves the RT problem self-consistently. Instead of the iterative procedure of 
estimating the dust temperature distribution (as described in Wolf et al. 1999), the concept 
presented by Bjorkman & Wood (2001), which is based on the immediate correction of the 
dust grain temperature after absorption of a photon package, has been used. Furthermore, 
the method described by Lucy (1999), which takes into account the absorption not only at 
the points of interaction of the photons with the dust grains but also in-between, has been 
applied. While the first method allows to simulate the RT for optical depths^ r ^ 10^, the 
second was used in order to increase the efficiency of the simulation for low optical depths. In 
addition to the test of MC3D described by Wolf et al. (1999) considering a single grain size, 
it has been successfully tested for the case of dust grain mixtures against the one-dimensional 
code of Chini et al. (1986). 

The RT is simulated at separate wavelengths within the wavelength range [Amin, Amax]- 
For this reason, the radiation energy of the emitting source(s) is partitioned into so-called 
weighted photons (Fischer, Henning, & Yorke 1999) each of which is characterized by its 
wavelength A and Stokes vector I — {I,Q,U,V)'^. In order to consider dust grain ensem- 



^Throughout the whole article the quantity r refers to the optical depth based on the extinction cross 
section of the dust grains at a wavelength of 550 nm. 
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bles instead of a single grain species, the RT concept as described by Wolf et al. (1999), 
Lucy (1999), and Bjorkman & Wood (2001) has to be extended in the following manner: 



1. The mean free path length / between the point of emission and the point of the first 
photon-dust interaction (and - subsequently - between two points of interaction) is 
given by 

rext,i = -ln(l-Z), (1) 



np 
i=l 



J=l k=l 



np 



Ah, Z = J]A/i. (2) 



i=l 



Here, Text,i is the optical depth along the path of the length /, fext,i is the corresponding 
numerical accumulated optical depth (Text,i — Text,! — * 0), tid is the number of chemically 
different dust grain species, ns is the number of particle radii being considered, and 
fi. is the spatial coordinate corresponding to the ith integration point along the path 
length /. The quantity p(r) is the number density at the spatial coordinate f, the 
quantity Cext is the extinction cross section, and np is the number of integration points 
along the path. Furthermore, Z is a random number uniformly distributed in the 
interval [0, 1]. 

2. According to the concept of immediate reemission described by Bjorkman & Wood (2001), 
at each point of interaction of a photon with a dust grain either scattering or absorp- 
tion occurs. The probability for a photon to undergo the one or the other interaction 
process with a grain with the chemical composition and the particle radius is 
given by 

jT — Pj,k{^ ^^j,k /o\ 

Z^m=l Z^p=l lPm,p{' J '^extm,p\ 

where 'x' stands either for absorption or scattering. In case of absorption, the imme- 
diate reemission occurs from the same dust grain (species). 

3. In Lucy's concept (1999), which considers the absorption of the electromagnetic radia- 
tion field not only at the end points of the photon path (points of interaction) but also 
in-between, simply the absorption due to all dust grain species (instead of only one) 
has to be taken into account: 

"^S /.point: 



'-a /•point2 

Tabs = J2J2 ^^bsi,k Pj,k{^ dr. (4) 

j=l k=l >^PO™ti 



The investigations presented in this paper are based on dust grain ensembles with dif- 
ferent grain size distributions but the same chemical composition. Since the formalism ( and 
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therefore at least the qualitative conclusions) are the same for ensembles consisting of grains 
with different chemical composition or size distribution or even both, this restriction is jus- 
tified. The optical properties of astronomical silicate (Draine & Lee 1984) for [A] a size 
distribution of small grains with radii a = (0.005 — 0.25) /im and [B] a size distribution of 
larger grains with radii a = (0.1 — 1) /im have been used in Sect. 3, 4.1, and 5. 

The grains are assumed to be spherical with a size distribution described by the widely 
applied power law n{a) oc a~^'^ (Mathis, Rumpl, & Nordsieck 1977). The optical properties 
such as the extinction and scattering cross sections as well as the scattering distribution 
function have been derived on the basis of Mie scattering, using the Mic scattering algorithm 
published by Bohren & Huffman (1983). The correct (anisotropic) scattering distribution 
function (for each dust grain size) has been considered in the RT process. 



In case of non-self-consistent RT in a dust grain mixture, the numerical effort (finally, 
the computing time and RAM requirement) may be substantially decreased by assuming 
weighted mean values of those parameters which describe the interaction of the electromag- 
netic field with the dust grains: the extinction, absorption, and scattering cross section, the 
Stokes parameters, and - as as function of the extinction and scattering cross section - the 
Albedo. 

The weight Wj{a), which represents the contribution of the jth component of the dust 
grain mixture, results from the abundance of this component in respect of its dust grain 
number density (assuming no chemically different dust species) and the size distribution of 
the respective material: 



where amin and Omax are the minimum and maximum grain radius of the size distribution. 
The Stokes parameters as well as the extinction, absorption, and scattering cross section 
(Cext, C'abs) C'sca) a^c additive. Therefore, the representative values in case of a dust grain 
mixture can be derived on the basis of their weighted contributions (see, e.g., Martin 1978, 



3. Mean dust grain parameters 




(5) 



Sole 1980): 
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and 



(Csca) = / Wj{a) Csc>,j{a) da, (8) 



{S) = J2 Wj{a) Sj{a) da, (9) 

where 5 is the Mueller matrix which is used to describe the modification of the photon's 
Stokes vector due to the interaction of a photon with the scattering/ absorbing medium (dust 
grains; see Bickel & Bailey 1985, Bohren & Huffman 1983). For the albedo A follows: 

Sj=i ^^^^^ ^^^^t^ 

Since this formalism (Eq. 5-10) covers both grain size and chemical distributions, the conclu- 
sions to be derived for a grain size distribution but a single chemical component is also valid 
for grain ensembles consisting of grains with different chemical compositions. The mean 
dust grain parameters used in the following sections have been derived on the assumption 
of 10^ grain radii equidistantly distributed in the grain size range of the small/large grain 
ensemble. 

In thermal equilibrium, the temperature of a dust grain (wth grain size or composition) 
can be estimated from the assumption of local energy conservation. If 

noo 

dEt^^dt Lf^dX (11) 

is the energy being absorbed by the uih. dust grain component during the time interval dt, 
then the local energy conservation can be written as 

= dt / dX^dt 47ra2 Qf^{a) 7rB^{T^)dX = d^f ^ (12) 

Jo Jo 

where dE'i^'^ is the total amount of energy being rc-cmittcd during the time interval dt, 
Lf^^ and U^*^ are the monochromatic absorbcd/rc-emitted luminosity at the wavelength A, 
Qx'ui'^) is absorption efficiency and B\{Ty) is the Planck function of the grain with the 
temperature T^. For two particles of the same composition but different radii (ai, 02) the 
temperature is consequently only the same if 

Qf\a,) = Qf\a,). (13) 



This is usually not the case. Thus, mean dust grain parameters may be used only in case of 
the simulation of dust scattering but not for the estimation of the dust grain temperature 
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distribution and resulting observables, such as the spectral energy distribution and images 
in the mid-infrared to millimeter wavelength range. 

In the following Sections (§4 and §5) it will be investigated under which conditions the 
approximation of the mean dust grain parameters provides nevertheless a good estimate 
for these observables. While the most crucial information can be derived on the basis of 
a one-dimensional model (§4), the study presented in §5 is aimed to reveal the influence 
of geometrical effects in case of two-dimensional disk-like structures. In each case, typical 
parameters for both the dust grain composition (see §2) and the circumstellar shells/disks 
are considered. 



4. Temperature Distribution 

4.1. Basic Model: Spherical Shell 

In order to reveal the differences between the results of the RT based on real dust grain 
size distributions and the approximation of mean dust grain parameters representing the 
weighted mean optical parameters of the grain size distribution, a simple-structured but in 
respect of the astrophysical applications very useful model has been used for the following 
investigations. It consists of a spherical shell (outer radius: 10^ AU; inner boundary: 0.2 AU) 
with a radial density profile described by a power law, p(r) oc r°, where a; is a negative, 
constant quantity. The dust grain size distribution is assumed to be constant throughout 
the whole shell. The only heating source is an embedded, isotropically radiating star in the 
center of the shell with an effective temperature and radius identical to that of the Sun. 

In order to study the influence of multiple scattering and (additional) heating due to 
the dust reemission, the optical depth of the shell r (as seen from the star) and the expo- 
nent «, describing the relative density gradient in the shell, have been considered as variable 
parameters of the model. Four values (-1.0, -1.5, -2.0, and -2.5) for the exponent a were 
chosen, covering a broad range of astrophysical objects ranging from circumstellar shells 
(see, e.g., Chini et al. 1986, Henning 1985, Yorke 1980) and disks (see, e.g., Men'shchikov & 
Henning 1997, Sonnhalter et al. 1995, Lopez et al. 1995) to the radial density profile being 
assumed for the dust density distribution in AGN models (see, e.g., Manske et al. 1998, 
Efstathiou & Rowan- Robinson 1995, Stenholm 1994). The optical depth r was varied from 
the optical thin case (t=0.1) to the case of intermediate optical depth (r = 1) and finally 
to r = 10 as the optically thick case. For optical depths below this interval no remarkable 
differences to the case r = 0.1 arc expected since multiple scattering and reemission have a 
negligible influence on the dust grain temperature which is determined only by the distance 
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from the star (attenuation of the stellar radiation field) and the absorption efficiency of the 
dust grains. For optical depths r 3> 10 preparatory studies have shown that the tempera- 
ture difference between the different dust grains becomes negligible, converging towards the 
temperature distribution obtained on the basis of the mean dust grain approximation. This 
finding is in agreement with results obtained by Kriigel & Walmsley (1984 - investigation of 
dust and gas temperatures in dense molecular cores), Wolfire & Churchwell (1987 - study 
of circumstellar shells around low- mass, pre- main-sequence stars), and Efstathiou & Rowan- 
Robinson (1994 - multi-grain dust cloud models of compact HII regions). Furthermore, 
only the inner region of the circumstellar shell shows a broad distribution of dust grain tem- 
peratures at a given distance from the star. In Fig. 1 and 2 the differences between the 
radial temperature distribution of the single dust grains T(r, a) and that resulting from the 
approximation of mean dust grain parameters T{r) is shown from the inner boundary of the 
shell to a distance of 1 AU from the star. The main results are: 

1. The temperature of the different grains spans a range of ~ 250 K (pa 25% of the 
corresponding mean temperature) around the temperature of the mean grains at the 
inner boundary of the shell. 

2. The temperature difference decreases 

• towards larger distances from the star, and 

• with increasing optical depth as soon as the shell becomes optically thick. 

3. The relativ density gradient, described by the exponent a is of minor importance. How- 
ever, an increase of the absolute amount of this exponent results in an increased redis- 
tribution of energy between the different dust grain components at the inner boundary 
of the shell (r > 1) and therefore in a decrease of the temperature dispersion. 

4. The difference {T{r, a) — T{r)) strongly depends on both the absorption efficiency, and 
thus on the dust grain size distribution, and the temperature of the heating source since 
the combination of these parameters determines the amount of energy being absorbed. 
This is clearly illustrated by the different signs of the temperature difference of the large 
and the small grain size distribution to the temperature of the mean grains shown in 
Fig. 1 and 2. 

The found complex, highly dispersed temperature structure in a dust grain mixture requires 
a large number (here: 2^) of considered grain radii to be correctly taken into account and 
cannot be adequately represented by a mean particle. This finding is of special interest for 
the simulation of chemical networks, since the existence of ice layers, the possibility of a 
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certain reaction on the grain surface, surface reaction rates, and the temperature of the sur- 
rounding gas phase depend on the dust grain temperature (see, e.g., Cazaux & Tielens 2002, 
Markwick et al. 2002, Charnley 2001). The dust phase in the close steUar environment is of 
particular importance, since the chemical evolution takes place here on its smallest timescale. 
Furthermore it was found, that the approximation of mean particle properties is found to be 
justified for dust configurations of high (r 3> 10) optical depth. 



4.2. HH30 circumstellar disk 

Although all basic characteristics of the temperature structure in a dust grain mixture 
can be studied on the basis of a one-dimensional model (Sect. 4.1), a disk-like structure 
will be considered in the following. The motivation for this is to investigate the possible 
influence of geometrical effects which have been minimized in case of the spherical shell. 
Furthermore, two-dimensional models are of particular importance since they are widely 
applied in simulations of circumstellar disks, debris disks around evolved stars, active galactic 
nuclei, and galaxies (see, e.g., Chiang & Goldreich 1997, 1999; see also the references in 
Sect. 4.1). 

In the following, a model of the nearly edge-on disk around the classical T-Tauri star 
HH30 is considered (see, e.g.. Burrows et al. 1996 for observational details). Following 
Cotera et al. (2001) and Wood et a. (2002), a flared geometry as described by Shakura & 
Sunyaev (1973) is adapted: 

(14) 




h{u) 



where cu is the radial coordinate in the disk midplane, po is the density at the stellar surface, 
and the scale height h increases with radius: 



ft = Ao^^j . (15) 

In order to reproduce the high density (and therefore temperature) gradient in the inner 
region of the disk, a very high resolution of the model both in radial and vertical direction 
has been applied. The smallest resolved structure at the inner radius of the disk has a 
hnear extent of about 29 % both in vertical and radial direction. As for the SED modeUing 
performed by Wood et al. (2002), the following values have been used: stellar radius i?* = 
1.2i?0, stellar effective temperature = 3500K, a = 2.25, P = 1.25, ho = 0.017i?*, inner 
radius of the disk Rq — 6i?*, and the outer radius of the disk amounts to 200 AU. In 
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contrast to the simulations presented by Cotera et al. (2001) and Wood et al. (2002), only 
one chemical dust grain component (astronomical silicate as in Sect. 4.1) has been chosen 
in order to simplify the simulation analysis. The size distribution is specified using a power 
law with exponential decay: 



(see, e.g., Kim et al. 1994). Following Cotera et al. (2001), a maximum grain size Omax — 
20//m has been chosen (other parameter values: Oc = 0.55//m, p = 3.0). A total number 
of 32 grain sizes with radii equidistantly distributed on a logarithmic scale between 0.05 //m 
and 20 yum are considered. The visual extinction of the disk amounts to Ay — 2800 at an 
inclination of 84° (sec Wood et al. 2002, Tabl. 2). The RT simulations have been performed 
on the basis of the model geometry shown in Fig. 6[B] in Wolf, Henning, & Stecklum 1999. 
While the star could be assumed to be point-like in the one-dimensional model (§4.1), its real 
extend has to be taken into account now. The radiation characteristic at each point of the 
stellar surface is described by a / oc cos(6') law (star without limb- darkening; see Wolf 2002 
for details). 

The resulting temperature distribution for two selected grain sizes (a = 0.05/im and 
a = 20/im) and their difference arc shown in Fig. 3. On the largest, ~10-100AU scale the 
temperature of the small grains is up to 40 K higher than that of the large grains in the 
optically thin atmosphere above/below the disk. In agreement with the results found in 
optically thick spherical shells, the temperature difference decreases towards the optically 
extremly thick midplane where it amounts to only a few Kelvin. 

However, the vertical temperature structure shows a more complex behaviour in the 
inner region of the disk, on a size scale of a few AU and smaller. Here, an inversion of 
the temperature difference in the upper layer of the optically thick disk occurs, i.e., the 
large grains are substantially warmer than the small grains. This remarkable effect can be 
explained by the more efficient heating of large grains by the mid-infrared to far-infrared 
reemission of the hot inner disk. This explanation is also supported by the comparison 
of the vertical temperature structure obtained at different distances from the star^: Fig. 4 
shows the difference between vertical temperature distribution of grains with increasing radii 



•^In order to obtain the spatially high-resolved temperature distributions shown in Fig. 4 and 5 on a PC 
with 1 GByte RAM, the RT has been simulated only in the inner disk region with a diameter of 8.8 AU, 
which is one to two orders of magnitude larger than the regions considered in the following discussions. The 
comparison with the results obtained on a temperature grid with a lower resolution (Fig. 3) showed that the 
resulting differences in the spatial temperature distribution is in the order of or smaller than the statistical 
noise of the results and therefore negligible. 




(16) 
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(a = 0.01 //m, 0.1 //m, 1.1 /xm, 9.0/xm) and a = 20 //m grains. Furthermore, Fig. 5 shows 
the temperature difference AT = T(a) — T{a = 20/im) as a function of the grain radius a 
and radial distance from the star at different distances from the midplane. While at a radial 
distance of 0.05 AU the temperature of the 20/im grains is about 130 K higher than that of 
0.01 /xm grains in the "inversion layer", the amount of this temperature difference drops to 
less than ~ 20K at a radial distance of 0.5 AU from the star (see Fig. 4). Consequently, the 
vertical extent and location of the minimum of the temperature inversion region depends 
on the particular grain size (and the vertical density distribution in the disk). For instance 
it was found that the minimum of the AT = T{a — 9.0/xm) — T{a — 20//m) distribution 
shown in Fig. 4 is shifted towards smaller distances from the midplane in all three considered 
cases (a)=0.05AU, 0.1 AU, and 0.5 AU). According to the explanation of the temperature 
inversion given above, a similar shift of the minima of the difference temperature distribution 
is expected in case of the other grain sizes (0.01 /xm to 1.1 /xm) as well. However, the spatial 
distribution of the simulated temperature distribution is too low to allow the verification of 
this assumption. 



5. Spectral Energy Distribution 

Based on the dust grain temperature distribution shown in Fig. 1, the spectral energy 
distribution (SED) for the different (spherical) shells have been calculated assuming ns = 2*s 
{is = 1, 2, . . . , 6) different grain sizes. The relative difference 

AS, = ^^^"l^ ~ (17) 

is shown in Fig. 6 and 7. Here, Sx{ns) represents the SED of the real grain size distribution, 
and Sx is the SED resulting from the approximation of mean dust grain parameters. 

As Fig. 6 illustrates, the results of the real grain size distribution converges towards 
the results based on the approximation of mean dust grain parameters. As it also shows, 
a minimum of 2^ to 2^ grain sizes have to be considered to achieve deviations of less then 
~ 10%. If too few single grain sizes are considered, the SED does not represent the observable 
SED of the real grain size distribution sufficiently but drastically overestimates 

1. The absorption of the stellar radiation (and therefore underestimates the visual to 
near-infrared SED), 

2. The depths of the absorption bands in general (as the silicate absorption feature at 
~ 9.7 /xm demonstrates), and 
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3. The dust reemission spectrum in the near-infrared to miUimeter wavelength range. 

While the different absorption behaviour influences the temperature structure and over- 
all SED in general, the influence of the different scattering behaviour depends on the optical 
depth. In the optically thin case, differences in the scattering behaviour result mainly in 
deviations of the short wavelength region of the stellar SED. At higher optical depth, the 
scattering behaviour has direct influence on the energy transfer in the dust envelope and 
thus on the temperature structure and the near- to far-infrared wavelength range as well. 
However, this effect is of importance only at intermediate optical depths (r 1 ... 10), since 
for even higher optical depths the temperature structure, at least in case of a one-dimensional 
dust configuration, tends not to depend on the description of the dust grain ensemble (real 
dust grain size distribution / approximation of mean dust grain parameters). 

6. Conclusions 

In simulations of the radiative transfer in the circumstellar environment of young stellar 
objects it has been widely established to use mean values for those parameters which describe 
the interaction of the electromagnetic field with the dust grains (see, e.g., modelling efforts 
by Calvet et al. 2002: circumstellar disk around the young low-mass star TW Hya; Cotera 
et al. 2001 and Wood et al. 2002: circumstellar disk around the classical T Tauri star HH 30; 
Fischer, Henning, & Yorke 1994, 1996: polarization maps of young stellar objects). A main 
reason for this is given by the fact that the consideration of a large number of single grain 
sizes and chemically different grain species results in a nearly linearly increasing amount of 
required computer memory in order to store the separate temperature distributions. Fur- 
thermore, the calculation of the spatial temperature distribution for all grains of different 
sizes and chemical composition requires significantly more computing power since (a) the 
heating by the primary sources (e.g., the star embedded in a circumstellar envelope) has to 
be performed independently for each grain species, and (b) the number of computing steps 
required to model the subsequent mutual heating of the different dust grain species due to 
dust re-emission scales even as ~ (n^ x ris)^, where rio is the number of chemically different 
components and ns is the number of separate dust grain radii considered. While these simu- 
lations are feasible in case of one-dimensional models (see, e.g., Chini, Kriigel, & Kreysa 1986 
and Efstathiou & Rowan-Robinson 1994: Dust emission from star forming region; Kriigel & 
Siebenmorgen 1994 and Siebenmorgen, Kriigel, & Zota 1999: Radiative transfer in galactic 
nulei) or simple-structured two-dimensional models (Men'shchikov & Henning 1997: Cir- 
cumstellar disks; Efstathiou & Rowan- Robinson 1994: Disks in active galactic nuclei), it is 
hardly possible to handle two- and three-dimensional models with high density gradients 
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and/or high optical depth which require high- resolution temperature grids. 

In this study, the difference between the results of RT simulations based (a) on a mean 
dust grain parameter approximation and (b) real dust grain size distributions have been 
investigated. Based on a one-dimensional density distribution it was found that the temper- 
ature structure of a real grain size distribution shows a very complex behaviour in the inner, 
hot region of the shell depending on (1) the grain size distribution, (2) the effective temper- 
ature of the embedded star and the optical depth and therefore on the density distribution. 
However, the relative difference between the SED based on a real dust grain size distribution 
on the one hand and the approximation of mean dust grain parameters on the other hand 
was found to be smaller than about ~10% if a minimum number of 2^ to 2^ grain sizes have 
been considered. 

As the temperature structure in a circumstellar disk - based on the model for HH 30 - 
shows, the geometry of the density distribution is a significant parameter for the resulting 
temperature differences between grains of different size, too. In the inner region of the disk 
with a diameter of a few AU a temperature inversion layer was found where the sign of the 
temperature difference of the largest and smallest grains is reverted. As this and the results 
obtained on the basis of the spherical shell (§4.1) show, the dust grain temperature structure 
is not sufficiently represented by dust grains with mean optical parameters. On the one 
hand, this is of tremendous importance for the simulation of chemical networks since the 
largest deviations from results based on the approximation of mean dust grain parameters 
have been found in the inner hot, dense region of the shell/disk where the chemical evolution 
takes place on its smallest timescale. On the other hand, the complex temperature struc- 
ture may significantly change the hydrostatic properties of the considered gas/dust density 
distribution itself. However, these questions have to be investigated in future studies in 
order to find out the influence on observable quantities such as SEDs, images, polarization 
maps, and visibilities. Furthermore, the influence on processes taking place in more evolved 
circumstellar disks, such as the dust settling and dust grain growth, have to be considered 
taking into account the temperature structure of real grain size distributions. 

This research was supported through the HST Grant GO 9160, and through the NASA 
grant NAGS- 11645. 1 wish to thank the referee E. Dwek who helped to improve the clarity 
of the presentation of the results. 
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Fig. 1. — Temperature difference of the dust grains T(a, r) — T{r) in the innermost region of 
the shell (inner boundary at 0.2 AU-1 AU). The quantity T{a,r) is the radial temperature 
distribution of the dust grains with the radius a; T{r) is the radial temperature distribution 
resulting from the approximation of the mean dust grain parameters. Small dust grain size 
distribution: a — (0.005 — 0.25) //m. The number of grain sizes being considered amounts 
to 2^ 




Fig. 2. — Temperature difference of the dust grains T{a, r) — T{r) in the innermost region 
of the shell (inner boundary at 0.2 AU-1 AU). Large dust grain size distribution: a = (0.1 — 
1) /im. Note that the temperature of the dust grains mainly increases with decreasing particle 
radius, which is in contrast to the behaviour of the small grain size distribution (see Fig. 1; 
see Sect. 4.1 for explanation). 
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Fig. 3. — Left and middle column: Color coded temperature structure of the HH30 circum- 
stellar disk model (left: a = 0.05/im dust component, middle: a = 20//m component). Right 
column: Temperature difference AT = T{a = 0.05 fim) — T{a = 20/im). The background 
color in the dust free regions above/below the disk represents the AT = level. Thus, re- 
gions with brighter colors than the background color represent areas where the large grains 
have a higher temperature than the small grains (see also Fig. 4). [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 4. — Temperature difference AT = T(a) — T(a = 20fiTa) as a function of tlie distance 
from the midplane z, obtained at radial distances (D=0.05AU, 0.1 AU, and 0.5 AU from tlie 
star. Considered particle sizes: a = 0.01 /im (solid line), 0.1 /im (dotted line), 1.1 /im (dashed 
line), and 9.0 ^m (dash-dotted line). 
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Fig. 5. — Temperature difference AT = T(a) — T(a = 20/um) as a function of the grain 
radius a obtained at radial distances (D=0.05AU, 0.1 AU, and 0.5 AU from the star. In each 
figure, the temperature difference AT is shown at 5 points with different distances from the 
midplane z which are located relative to the coordinate z = at which the minimum of 
the temperature difference for a = 0.01 fim, 0.1 /xm (dotted line), and 1.1 fim was found in 
Fig. 4: z = O.Szm (dotted line), z = 0.75z^ (short-dashed line), z = z-^i (dot-dashed line), 
z — 1.25ztn (dot-dot-dot-dashed line), and z — l.bz^a (long-dashed line). 
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Fig. 6. — Relative difference ASx of the SEDs resulting from a real dust grain mixture based 
on 2*^ (is = 1, 2, . . . , 6) dust grain sizes and SEDs based on the approximation of mean dust 
grain parameters. See the upper left Figure for the explanation of the different drawing 
style. A close-up view on the difference spectrum ASx in the case of 2^ grain sizes is shown 
in Fig. 7. 
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Fig. 7. — Relative differences AjS'a between tlie SEDs based on 2^ dust grain sizes and 
tlie SED based on tlie approximation of mean dust grain parameters. Tlie influence of the 
exponent a (see §4.1) is clearly of minor importance. 



